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Abstract 

In this manuscript we give thought to the aftermath on the stable probability density 
function when standard multiplicative cascades are generalised cascades based on the 
g-product of Borges that emerged in the context of non-extensive statistical mechanics. 

1 Introduction 

In the twenty years that have elapsed since the publication of the non-additive entropy 
S q , also fairly known as Tsallis entropy [I], many applications and connections to natural 
and man-mind phenomena have been established [5J. One of the most exciting which 
that have emerged within the non-extensive scope is the definition of a whole new set of 
mathematical operations/functions that goes from the generalised algebra independently 
defined by Borges [3] and Nivanen et al. [I] and the integro-differential operators by Borges 
to the g-trigonometric functions [5]. Besides its inherent beauty, these generalisations 
have found its own field of applicability. Namely, the g-product plays a primary role 
in the definition of the g- Fourier transform [6], thus in g-Central Limit Theorem [7], 
whereas the generalised trigonometric functions have been quite successful in describing 
the critical behaviour of a class of composed materials known as manganites |8j. In 
this article, we inquire into the possible applications of the g-product in the generation 
of random variables and its consequence on the definition of a new class of probability 
density functions. 

2 Preliminaries: the g-product 

The g-product, <S> g , has been introduced with the purpose to find a functional form that 
is able to generalise in a non-extensive way the mathematical identity, 

exp [In x + In y] = x x y, (x, y > 0) , (1) 

so that the equality, 

x <S) q y = exp q [ln q x + \n q y] , (2) 
holds. The representations ln g (.) and exp g (.) correspond to the g-logarithm [9], 

\n q (x) = X \ — , (x>0,gGK), (3) 
1 — g 
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and its inverse, the g-exponential, 

exp g (x) = [1 + (1 -q) x]^ , (x,qeft), (4) 

respectively (exp g (x) = if 1 + (1 — q) x < 0). For q — > 1, the equation ([2]) recovers 
the usual property, 

In (x x y) = In x + In y 

(x, y > 0), with x x y = x <8>i y. Its inverse operation, the g-division, x g verifies the 
following equality (x ® q y) <2) q y = x. 

Bearing in mind that the g-exponential is a non-negative function, the g-product must 
be restricted to the values of x and y that respect the condition, 

\x\ x ~ q + M 1-9 - 1 > 0. (5) 

Moreover, we can extend the domain of the g-product to negative values of x and y writing 
it as, 

x® q y= sign (xy) exp g [In, |x| + ln 9 |y|] . (6) 
Regarding some key properties of the g-product we mention: 

1. x <8>i y = x y; 

2. x ® q y = y ® q x; 

3. (x ® q y) ® q z = x <g> 9 (y ® 9 z) = [x l - q + y l ~ q - 2]^; 

4. (x (g) q 1) = x; 

5. hiq [x ®g 2/] = hiq x + ln g y; 

6. ln g (xy) = ln q (x) + \n q (y) + (1 - q) hx q (x) In, (y); 

7. (X® g ?/) _1 = X" 1 ®2-gZ/ _1 ; 

{0 if (g > 1 and x > 0) or if (g < 1 and < x < 1) 

(x 1-9 — l) 1-9 otherwise 

For particular values of g, e.g., g = 1/2, the g-product provides nonnegative values at 
points for which the inequality |x| 1-9 + \y\^ q — 1 < is verified. According to the cut-off 
of the g-exponential, a value of zero for x ® q y is set down in these cases. Restraining our 
analysis of the Eq. (jSJ) to the sub-space x,y > 0, we can observe that for q — > — oo the 
region {0<x<l,0<y<l} is not defined. As the value of q increases, the forbidden 
region decreases its area, and when q = 0, we have the limiting line given by x + y = 1, for 
which x®o2/ — 0. Only for q = 1, the entire set of x and y real values of has a defined value 
for the g-product. For g > 1, the condition (JSJ) implies a region, |x| 1-9 + = 1 for 

which the g-product diverges. This undefined region augments its area as q goes to infinity. 
When q = oo, the g-product is only defined in {x > 0, < y < 1} U {0 < x < 1, y > 1}. 
Illustrative plots are presented in Fig. (1) of Ref. [TO] . 
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3 Multiplicative processes as generators of distribu- 
tions 



Multiplicative processes, particularly stochastic multiplicative processes, have been the 
source of plenty of models applied in several fields of science and knowledge. In this con- 
text, we can name the study of fluid turbulence [UJ, fractals [12], finance [JjJj, linguistics 
[H] , etc. Specifically, multiplicative processes play a very important role on the emergence 
of the log-Normal distribution as a natural and ubiquitous distribution. In simple terms, 
the log-Normal distribution is the distribution of a random variable whose logarithm is 
associated with a Normal distribution [To] . 



1 
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(7) 



With regard to the dynamical origins of the log-Normal distribution, several processes 
have been thought up to generate it. In this work we highlight the two most famous 
- the law of proportionate effect [16] . the theory of breakage [TT] or from Langevin-like 
processes [IB] . We shall now give a brief view of the former; Let us consider a variable Z 
obtained from a multiplicative random process, 

N 

Z = ]JCi, (8) 

i=i 

where Q are nonnegative microscopic variables associated with a distribution /' \ £j . If 
we consider the following transform of variables Z = In Z, then we have, 

N 



i=l 



with £ = ln^. Assume now ( as a variable associated with a distribution /(C) with 
average [i and variance a 2 . Then, Z converges to the Gaussian distribution in the limit of 
N going to infinity as entailed by the Central Limit Theorem [19J. Explicitly, considering 
that the variables ( are independently and identically distributed, the Fourier Transform 
of p (Z') is given by, 



T\p{Z')\ (k) 



e lk -f(C) dC, 



N 



(9) 



where Z' = N l Z . For all N, the integrand can be expanded as, 

~l N 

n! N 



T\p{Z')\ (k) 
T\p{Z')\ (k) = exp ijVln 



oo 

E- 

?1=0 



(10) 



1 + ikf -hk 2 ^- 1 



AT2 



+ 0(N~ 
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expanding the logarithm, 



T [P (Z')\ (fc) « exp 



11) 



Applying the inverse Fourier Transform, and reverting the Z' change of variables we finally 
obtain, 



p(Z) 



1 



V2nNa 



exp 



2a 2 N 



(12) 



We can define the attracting distribution in terms of the original multiplicative random 
process, yielding the log-Normal distribution 



1 



V2 n No Z 



exp 



(inZ-Nfif 
2 a 2 N 



(13) 



Although this distribution with two parameters, /i and a, is able to appropriately 
describe a large variety of data sets, there are cases for which the log-Normal distribution 
fails statistical testing [JS]. In some of these cases, such a failure has been overcome by 
introducing different statistical distributions (e.g., Weibull distributions) or changing the 
2-parameter log-Normal distribution by a 3-parameter log-Normal distribution, 



p [Xj 



2 7TO" (X — 



exp 



(In [x-9]- fjf 
2o~ 2 



(14) 



In the sequel of this work we present an alternative procedure to generalise the Eq. ((7j). 
The motivation for this proposal comes from changing the N products in Eq. (jSJ) by N 
q-pro ducts, 

N ^ 

] [ Ci = Ci® 9 C 2 ® g ---® 9 Cjv- (15) 



i=l 



Applying the g-logarithm we have a sum of iV terms. If every term is identically and 
independently distributed, then for variables Q = ln q Q with finite variables we have a 
Gaussian has stable distribution, i.e., a Gaussian distribution in the g-logarithm variable. 
From this scenario we can obtain our q-log Normal probability density function, 



p [x) 

with the normalisation, 



1 



Z q xi 



exp 



(ln 9 x — p,)' 
2^2 



(x > 0) , 



(16) 



ferfc 



|erfc 



V2a V l ~1 
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V2a 



a if q < 1 
a if q > 1. 



(17) 
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In the limit of q equal to 1, ln g _»i x = lnx and Z q ^\ = v2 na and the usual log-Normal is 
recovered thereof (erfc stands for complementary error function). Typical plots for cases 
with g = |,g = l,g = | are depicted in Fig. [TJ 
The raw statistical moments, 

POO 

(x n ) = / x n p(x) dx, (18) 
Jo 



can be analytically computed for q < 1 giving 

T [u] exp 



7 

8/i 



with 



/3 



a/^ttct (1 - q) erfc 
1 + M (1 



\2 ! 
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(19) 
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2a 2 (i - g y 

where D_ a [z] is the "parabolic cylinder junction [21]. For q > 1, the raw moments are given 
by an expression quite similar to the Eq. (fl9l) with the argument of the erfc replaced by 

+ Z 1 ) • However, the fmiteness of the raw moments is not guaranteed for every 

q > 1 for two very related reasons. First, according to the definition of D^ p [z], v must 

be greater than 0. Second, the core of the probability density function, exp 

does not vanish in the limit of x going to infinity oo, 



(In, x-fj,) 



2 a 2 



lim exp 

This means that the limit p (x — > oo 



(ln 9 x-n)' 



2 a 2 









" I 2 ' 


= exp 


2 _ 



(21) 



is introduced by the normalisation factor x 



which comes from redefining the Gaussian of variables, 

V = hi g x, 



(22) 



as a distribution of variables x. Because of that, if the moment surpasses the value of q, 
then the integral (ITS]) diverges. 



4 Examples of cascade generators 

In this section, we discuss the upshot of two simple cases in which the dynamical process 
described in the previous section is applied. We are going to verify that the value of q 
influences the nature of the attractor in probability space. 

4.1 Compact distribution [0,6] 

Let us consider a compact distribution for indentically and independently distributed 
variables x within the interval and b. Following what we have described in the preceding 
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section, we can transform our generalised multiplicative process into a simple additive 
process of y^ variables which are now distributed in conformity with the distribution, 

p'{y) = \[l + {l-q)v]^ , (23) 

with y defined between -^j and if q < 1, whereas y ranges over the interval between 

— oo and ^j^" 1 when q > 1. Some curves for the special case 6 = 2 are plotted in Fig. [21 
If we look at the variance of this independent variable, 

°l = (y 2 ) - M 2 . ( 24 ) 

which is the moment whose finitude plays the leading role in the Central Limit Theory, 

3 
2 - 



we verify that for q > ~, we obtain a divergent value, 



2 b 2 - 2 * , s 

ol = 25 

y (3-2g)(2-g) 2 

Hence, if q < |, we can apply the Lyapunov's central Limit theorem and our attractor 
in the probability space is the Gaussian distribution. On the other hand, if q > |, 
the Levy-Gnedenko's version of the central limit theorem [22] asserts that the attracting 
distribution is a Levy distribution with a tail exponent, 

a = . (26) 

q-1 V ' 

Furthermore, it is simple to verify that the interval (§,oo) of q values maps onto the 
interval (0, 2) of a values, which is precisely the interval of validity of the Levy class of 
distributions that is defined by its Fourier Transform, 

T [L a (Y)] (jfe) = exp [-a \k\ a ] . (27) 

In Fig. [3] we depict some sets generated by this process for different values of q. 



4.2 q-\og Normal distribution 

In this example, we consider the case of generalised multiplicative processes in which the 
variables follow a (/-log Normal distribution. In agreement with what we have referred to 
in Sec. El the outcome strongly depends on the value of q. Consequently, in the associated 
x space, if we apply the generalised process to iV variables y = \n q x (x G [0, oo)) which 
follow a Gaussian-like functional form with average /i and finite standard deviation cr, 
i.e., V g< i or q > 3 in Eq. lfTB"]) . the resulting distribution in the limit of iV going to infinity 
corresponds to the probability density function (1161) with /i — > iV// and cr 2 — > iV a 2 . In 
respect of the conditions of q we have just mentioned here above, the g-log normal can be 
seen as an asymptotic attractor, a stable attractor for q — 1, and an unstable distribution 
for the remaining cases with the resulting attracting distribution being computed by 
applying the convolution operation. 

2 Strictly speaking, we cannot use the term Gaussian distribution because it is not defined in the 
interval (—00,00). The limitations in the domain do affect the Fourier transform and thus the result of 
the convolution of the probability density function. 
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5 Final remarks 



In this manuscript we have introduced a modification in the multiplicative process that 
has enabled us to present a modification on the log-Normal distribution as well as other 
distributions with slow decay. This distribution is controlled by an extra-parameter, q, 
when it is compared with the regular 2-parameter log-Normal distribution, that can be 
dynamically related to a change in the multiplicative random process. Besides, it provides 
interesting mechanisms of on-off dynamics. 

Regarding further applications, it is known that the standard log-normal distribution 
is unfitted for several data sets. This 3-parameter log-Normal probability function is 
expected to provide a better approach to these data [23J. 

SMDQ thanks his colleagues at Unilever for discussions and financial support of the Marie 
Curie Fellowship programme (European Union). 
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Figure 1: Plots of the Eq. ( fl6i) vs x for q = | (dotted line), q = 1 (full line) and q = | 
(dashed line) in linear-linear scale (upper), log-linear (centre), log- log (lower). 
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Figure 3: Sets of random variables generated from the process (1151) with N = 100 and 
q — — ~ (green), (red), ~ (blue), 1 (black), | (magenta) in linear (upper panel) and log 
scales (lower panel). The generating variable is uniformly distributed within the interval 
[0, 1] as is the same for all of the cases that we present. As visible, the value of q deeply 
affects the values of = Z. 
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